Remarks on the n = i disclination line in 
Landau-de Gennes theory of nematic liquid 

crystals 



H. Arodz and R. Peika 
Institute of Physics, Jagellonian University, 
Reymonta 4, 30-059 Cracow, Poland 



Abstract 

Using Landau-de Gennes effective theory for nematic liquid crystals 
we analyse structure of a rectilinear, n = 5, smooth disclination line 
in the case of equal elastic constants. We find that at certain tempera- 
ture there is an exact mathematical correspondence with a rectilinear 
vortex in superfluid 4 He. With a help of polynomial approximation 
difference of free energies of smooth and singular disclination lines is es- 
timated. It turns out that the smooth disclination line is energetically 
preferred only if temperature is low enough. At higher temperatures 
a disordered core should be expected. 

PACS numbers: 61.30.Jf, 11.27.+d 
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1 Introduction 



Disclination lines in nematic liquid crystals are among the most popular 
topological defects encountered in condensed matter physics ||], ||, 0. 
In spite of that, theoretical description of them still poses rather interesting 
problems. In the director formalism for uniaxial nematic liquid crystals there 
is the question of structure of a singular core of the disclination line. As it 
has been pointed out in papers ||, ||, [0], it is also possible to have smooth 
disclination lines, provided one allows for biaxiality induced by torques which 
are due to topologically nontrivial boundary conditions. This type of discli- 
nation lines has been investigated in the framework of Landau-de Gennes 
effective theory for nematic liquid crystals, mainly with the help of numer- 
ical methods. Analytic description, exact or approximate, is still missing 
except for special cases. Yet another theoretical description of the disclina- 
tion lines can be found in [§], where Ericksen effective theory is used, or in || 
where the biaxiality is not introduced. There are also less fundamental but 
more difficult problems, like theoretical description of evolution of a curved 
disclination line, see, e.g., ||TD|| . Actually, the present work is the first step in 
our attempt to solve the latter problem. 

We consider the rectilinear n = | disclination line within the framework 
of Landau-de Gennes effective theory. Review and discussion of this theory 
is given in, e.g., |lTJ. In the present paper we assume that the constant L 2 



is equal to zero, what corresponds to equal elastic constants in the director 
formalism. The case of L 2 ^ is much more complicated |12|] . First, we 
reformulate the theoretical description of the smooth disclination line found 
in [0], in particular in order to prepare a convenient starting point for the 
analysis of the curved smooth disclination line. Next, we present certain new 
results. We find that there is a special case in which the transverse profile 
of the smooth disclination line is given by one nontrivial function, instead 
of two in the generic case. We provide an approximate analytic description 
of the disclination line by means of the polynomial approximation. Finally, 
comparing free energies of the smooth and singular disclination lines we notice 
that as we vary parameters of the model each of the two can be energetically 
preferred. 

The contents of our paper is as follows. In Section 2 we recall formulas 
from Landau-de Gennes theory which we need in our considerations. Section 
3 is devoted to the analytic description of the smooth rectilinear disclina- 
tion line. This Section contains our new results mentioned in the preceding 
paragraph. In Section 4 we have collected several remarks. 
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2 Landau-de Gennes theory 



The order parameter for the nematic liquid crystal in Landau-de Gennes 
theory has the form of a symmetric, traceless real tensor Q = [Qij], where 
i, j = 1, 2, 3. The corresponding free energy density is given by the following 
formulas: 

1 1 

T = -LidiQ jk diQ jk + -LidiQikdjQjk + V(Q), (1) 

where 

V(Q) = -^Tr(Q 2 ) - b -Tr(Q 3 ) + ^(Tr(Q 2 )) 2 . (2) 

The constants a,b,c,Li,L 2 characterise the liquid crystalline material, see, 
e.g., |]T3| . In the following we assume that all the constants in (1), (2) are 
positive. Then, the ground states obtained by minimizing T are uniaxial, in 
accordance with what has been experimentally found for most of nematics. It 
is clear that in the ground state Q is constant. One can show by minimizing 
V(Q) that all ground states Q g can be obtained by uniform rotations or by 
reflections from 

/ 2 \ 
Qo = Vo -1 , (3) 
\0 -I j 

where 



b + Vb 2 + 2Aac 
^ = Uc • (4) 



That is, 



Q g = OQ O T , (5) 

where O G 0(3) is a constant orthogonal matrix. The corresponding value 
of V is equal to 

^min = ->7o(l + /?)> (6) 

br]o 
3a 



3 

mm ~ ~ 2 
where 

P = 



By definition, in the uniaxial case the matrix Q has two identical eigen- 
values Ai = A2 = A. Then one can define the director field n by writing Q in 
the form 

Q = 3X (h-nn T ), (7) 

where I denotes the 3 by 3 unit matrix, the vector n has unit length, and 
nn T is a matrix (dyad). It is clear from formula (7) that as the eigenvectors 
corresponding to the degenerate eigenvalue A we can take any pair of linearly 
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independent vectors orthogonal to n. The third eigenvector is just n and the 
corresponding eigenvalue is equal to — 2A. For the particular Qq given by 
formula (3) the director field can be taken in the form 



n = , A = —770, 




and 



n g = On (8) 



corresponds to Q g given by formula (5). 

The true value of the Landau-de Gennes theory is of course related to the 
fact that it also describes states which are not the ground states. Then Q can 
be space and time dependent, Q = Q(f,t). For example, one can consider 
time evolution of inhomogeneous perturbations of the ground states. Another 
class of states contains topological defects. These states are often discussed 
in London's approximation |L4| which consists in restricting all configurations 
to the set of minima of V(Q), often called the vacuum manifold. Thus, in 
this approximation Q has the form (7) where now n = n(f, t) while A = — 770 
remains constant. In London's approximation free energies of various states 
can differ only by the derivative terms in (1). Formulas (1) and (7) give the 
well-known Oseen-Zocher-Frank free energy 

T = \K n (divn) 2 + \K 22 (n ■ rotn) 2 + \K^(n x rotn) 2 

+\K 22 d\\[(n • V)n - ndivn] + K mm , (9) 

where 

K u = K 33 = 18r ] 2 {L 1 + ^L 2 ), K 22 = ltofiL 1 , (10) 

and the constant V m [ n is given by formula (6). From formulas (10) we see 
that L 2 = corresponds to K\\ = K 22 = K 33 . 

Unfortunately, in the case of disclination lines London's approximation 
is too restrictive. The point is that it allows only for disclination lines with 
a singular core: due to the topologically nontrivial boundary conditions Q 
necessarily leaves the vacuum manifold at certain points in the space. On 
the other hand, in the full framework of Landau-de Gennes theory also a 
smooth disclination line, which does not contain any singular core, is possible. 
Which of the two types of disclination lines is expected to occur in a concrete 
nematic material at a given temperature can be found out by checking the 
corresponding values of the free energy. 
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3 The smooth n = \ disclination line 
3.1 The axially symmetric Ansatz 

Topological charge pertinent to the n = \ disclination line is related to the 
fundamental group of connected components of the vacuum manifold, that 
is to iri(SO(3)/H), where H C SO (3) is the stability group of Q . It is a 
well-known fact that iti(SO(3)/H) = Z 2 , |T]J. Far away from the rectilinear 
disclination line, that is at the spatial infinity in directions perpendicular to 
the line, the order parameter Q lies in the vacuum manifold (5). The matrices 
O in formula (5) taken along a large circle around the disclination line should 
form a continuous path in the SO (3) group such that its projection on the 
vacuum manifold SO(3)/H gives a noncontractible loop corresponding to the 
nontrivial element of tti (SO (3) /H) . Let us assume that the disclination line 
is perpendicular to the (x, y) plane. Then we may take 



(9(0) 



/ cos | — sin I 

sin § cosf | . (11) 



V 



where (p, 4>) are the polar coordinates in the (x, y) plane. 

The structure of the disclination line is determined from the requirement 
that the total free energy 

' tfxT (12) 



has a minimum within the chosen topological class of the order parameter 
Q. The necessary condition for that has the form 

SF 

+ Xij - Xji + \5ij = 0, (13) 



SQii(x 



where A^-, A are Lagrange multipliers corresponding to the conditions Qij — 
Qji = 0, Qkk = 0, respectively. In the discussed case of equal elastic 
constants one may expect that the minimal free energy is obtained for an 
axially symmetric configuration. Therefore, we look for a smooth solution of 
Eq.(13) assuming the following cylindrically symmetric Ansatz 



( S(p) + 3R{p) 

S(p)-3R(p) |0 T (0), (14) 

V o o 




where f]o is given by formula (4), 0(4>) has the form (11), and S(p), R{p) are 
unknown functions of the polar radius p. Ansatz (14) can be written also in 
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the form 



^ 2 





( 1 







1 COS0 


sin 





S{p) 





1 


i) 


+ 3R{p) sin 


— COS 







v° 







V o 









(15) 



This Ansatz is equivalent to the one considered in fl7|, but the form (14) is 
more transparent from the homotopy group viewpoint. 

The appropriate boundary conditions for the functions S(p),R(p) have 
the form 

R(0)=0, R(oo) = l, (16) 

S(0)=w o , 5(oo) = l, (17) 

where Wq is a constant. They follow from continuity of Q at p = 0, and from 
the assumption that in the limit p — > oo Q has the form OQqO t with Q 
given by formula (3). It is clear from the Ansatz (14) and from the conditions 
(16), (17) that Q is not restricted to the vacuum manifold. 
The free energy density (1) is now equal to 
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S' 2 + 3R' 2 + jzR 2 , 



%S 2 - 2R 2 



-6f3S{R 2 - ±S 2 ) + Ml + 3(3){S 2 + 3R 



2\2 



(18) 



where ' stands for d/ds. For convenience we have introduced the dimension- 
p/fo, where ^ : 



less variable s = p/£o, where = y2Li/3a. 

For Q of the form (15) Eq.(13) reduces to the following two equations 
(recall that L 2 = 0) 



S" + -S' + -S --(1 + 3(3)S(S 2 
s 3 6 



3R 2 



(5{S 2 - 3R 2 ) = 0, 



1 



R" + -R' 

s 



\r + Ir 



1 



0. 



(19) 
(20) 



, . '1 + 3(3)R(S 2 + 3R 2 ) + 2(3SR 
s* 3 6 

Equations (19), (20) are quite complicated. One can notice one simple 
interesting solution of them besides the trivial R—S — 0, namely R — 0,S — 
—2. This solution describes the ground state configuration equivalent to Q . 
Another simple solution, namely R — 0, S — 2/(1 + 3/3), is not interesting 
because it corresponds to a local maximum of V(Q). 

Less trivial solutions can be studied with the help of numerical methods. 
In Figs.l, 2 we present two examples of numerical solutions of Eqs.(19), (20) 
with the boundary conditions (16), (17). The solutions are represented by 
the continuous lines. They were generated with the help of Maple©. The 
dotted lines in these figures represent approximate forms of the functions 
R, S obtained in Subsection 3.3 below. 
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Fig.l. The functions R, S for j3 = 0.1. The numerical solutions are repre- 
sented by the continuous lines, while the dotted lines represent the approxi- 
mate forms of R, S constructed in Subsection 3.3. For the numerical solution 
w = 1.35710(4). 
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3.2 The connection with Ginzburg-Pitaevskii equation 

The boundary conditions (16) imply that R(s) can not be constant, while 
for S(s) this is not excluded. Actually, Figs. 1 and 2 suggest that for certain 
value of /3 between 0.1 and 1 the function S is constant. If we assume that S = 
const, then Eq.(19) splits into two independent algebraic equations (because 
for R 2 we can substitute, e.g., and 1): 

^(1 + 3/3)5-/3 = 0, ^(l + 3[3)S 2 + [3S-^ = 0. 
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They imply that 

P=\, S=l, (21) 

and consequently 

b 2 

ac=-. 

Now the remaining equation (20) has the form 

R" + -R' -\r + R(l - R 2 ) = 0. (22) 
s s 2 

Precisely this equation appears also in the theory of superfluid 4 He, namely 
Ginzburg-Pitaevskii equation for a rectilinear vortex with the unit winding 
number can be written exactly in the form (22), ||16|| . Certain similarity 
between the smooth disclination line and the superfluid vortex has already 
been noticed in JDfl, where purely uniaxial nematic liquid crystal with Q 
restricted to the form (7) is considered. Now we see that in full Landau-de 
Gennes theory the pertinent equations just coincide. 

In the superfluid 4 ife case the vacuum manifold can be identified with 
the Abelian U(l) group, hence it is quite different from SO(3)/H. In spite of 
that, the coincidence of equations suggests that dynamical properties of the 
supefluid vortex and of the disclination line can be quite similar. 

The material constant a linearly depends on temperature T, namely 

a = a (^* — T). 
The nematic phase exists in certain finite temperature range, 

T„ > T > T m . 

The parameter (3 monotonically decreases from oo to (3 m when T decreases 
from T* to T m . 

For MBBA nematic liquid crystal the material constants have the values 
a = 42 x 10 3 Jm- 3 K-\ b = 64 x 10 4 Jm" 3 ,c = 35 x 10 4 Jm" 3 , and T* = 



320K,T m = 294K, see, e.g., [13|]. It follows from the definition of (3 given 
below formula (6) that (3 m ~ 0.174. The special value (3 = 1/3 corresponds 
to the temperature Ti/ 3 ~ 311.fr. This estimate of Tu s for MBBA has a 
blemish: in reality the elastic constant L 2 ^ 0. Therefore, the obtained 
value should be taken with a grain of salt. 
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3.3 The polynomial approximation 

Even for relatively simple Eq.(22) the relevant exact analytic solution is not 
known. Below we present an approximate solution of Eqs.(19), (20) obtained 
with the help of so called polynomial approximation |fnj| . This approximation 
has turned out to be very useful in several cases. It was compared with good 
results with purely numerical solutions in [18 , where a vortex in Abelian 
Higgs model was considered. 

We already know from Figs.l, 2 the general shape of the solutions R, S. 
The pictures suggest the following approximate description of these functions. 
For small s we approximate R and S by the low order polynomials, 

i?< = r lS , (23) 

S< = w + w 2 s 2 , (24) 

which can be regarded as truncated series solutions of Eqs.(19), (20). Equa- 
tion (19) gives the following recursive relation 

91 

(25) 



1 

w 2 = -wo 



hi + 3f3)w 2 + f3w - I 
.o 3 



Here ri and wo are unknown parameters as yet. On the other hand, when s 
is large we use the following approximate solutions of (19), (20), 

_ 1 + 15/3 1 1-3(3 1 

K> ~ 1 12(3(2 + 3(3)^ b> ~ + 4(3(2 + 3(3)^- {2b) 

The exact solutions are expected to be smooth for all s. Therefore, we 
postulate that i?< smoothly matches _R> at certain point si. To determine the 
two unknown parameters r±, S\ we may impose the two matching conditions 

i?<( Sl ) = R^s,), R'^st) = R'^sx). (27) 

They give 

2 ( 1 + 15/3 V /2 

ri = s? s ' = {m2TW)) ■ (28) 

Similarly, for S we impose the following two matching conditions at s = s 2 

5<( S2 )=5>( S2 ), S' < (s 2 ) = S' > (s 2 ), (29) 

from which we determine Wq and s 2 . Specifically, we obtain the following 
relation between wq and s 2 



1-3/3 1 
2/3(2 + 3(3) 7 2 



w o = 1 + , om 3» ( 3 °) 
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and the algebraic cubic equation for wq 



(1 - 9(3 2 )w 3 + 54/3(1 + (3)wl - 4(1 + 21/3 + 36/3 2 )w + 24/3(2 + 3/3) = 0. (31) 

This last equation can easily be solved numerically for each given value of (3. 
It is clear from formula (30) that if0</3<l/3 the relevant solution wq of 
Eq.(31) should be greater than 1, while w < 1 for /3 > 1/3, and of course 
wo — 1 for /3 = 1/3. We have obtained the following approximate formulas 
for S2 and Wo: 
when /3 ^> 1 

* 2 -^(l + ^ + 0(/T 2 )), - 0.635 + ^ + (9(/3- 2 ), (32) 
when /3 — > 

s 2 « ^= (l + ^ + 0(/3 2 )), W0 = 2-12/3 + O(/3 2 ). (33) 

In this way we have obtained the approximate forms of R, S in the full 
range of the s variable. They are presented in Figs.l, 2 (the dotted lines) for 
/3 = 0.1 and /3 = 1.0. In particular, the approximate values of wo following 
from Eq.(31) are equal to 1.36360(9) and 0.79234(4), correspondingly. They 
are in quite good agreement with the values for the numerical solutions given 
in the figure captions. 

The approximate solutions R, S are continuous together with their first 
derivatives (the C 1 class) for all s, they converge to the exact solutions at 
small and large s, and they obey the boundary conditions (16), (17). There- 
fore, these approximate forms of R, S give a smooth disclination line which 
belongs to the right topological class. Its free energy (in a finite volume) is 
a little bit larger than the one for the exact solutions of Eqs.(19), (20). 

The Ansatz (14) and our approximate solution imply that at the center 
of the smooth disclination line 

Tr(Q 2 )( S = 0) = i«, 
while in the ground state (3) 

Tr(Ql) = Q V 2 . 

We see that degree of ordering (i.e. Tr(Q 2 )) inside the smooth disclination 
line is smaller than in the ground state. Only in the limit (3 — > 0, when 
w — > 2, the degree of ordering approaches its ground state value. 
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Let us end this Subsection with two remarks about the polynomial ap- 
proximation. First, there is a caveat in constructing the approximate forms 
of R and S. If we include more terms in the expansions (23-26), then it turns 
out that always R> > 0, while the truncated series R' K has coefficients with 
alternating signs. Therefore, it may happen that the matching conditions 
(27) can not be satisfied if _R< contains a wrong number of terms: at the 
matching point s\ the highest order term in the polynomial R' K can domi- 
nate and if the sign is minus the matching is not possible at all. In the case 
at hand it turns out that we may take either first or fifth order polynomials, 
but not of the third order one. Analogously, one has to be cautious when 
choosing a polynomial for S < . 

Second, including more and more terms in the polynomials R <7 S< in 
general does not seem to be the most efficient way to improve the approxi- 
mation. The point is that the infinite series expansions for R, S in s at small 
s, and in s -1 at large s, can be slowly convergent, or even they can have 
a finite and too small radius of convergence. We think that it is better to 
take the polynomials of the lowest possible order, to construct the approxi- 
mate solution, and then to compute corrections 5R } 5S to it by linearising 
equations (19), (20) around the approximate solution. The reason is that 
the approximate solution constructed above already takes care of the non- 
trivial, nonlinear structure of the disclination line, including the boundary 
conditions (16), (17), and therefore we expect that SR, 5S will be small. The 
solutions of the resulting linear equations for 5R, 5S are accessible by many 
techniques. 



3.4 Free energy of the smooth disclination line 

The approximate formulas for the functions R and S can be used to estimate 
the free energy per unit length of the smooth disclination line. The free en- 
ergy density decreases at large s rather slowly: the polynomial approximation 
inserted in formula (18) gives for s — > oo 

T ^y . (34) 
J *min ~ g 2' v ) 

where V m ^ n is given by formula (6). This asymptotic behaviour is the same 
as for the singular disclination line. General shape of the free energy density 
can be seen from Fig. 3, where we have plotted /(s) = 8(jF(s) — V m ^ n ) / (9ar]Q) 
for = 0.1. 
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2 4 6 8 10 12 14 "?6 

Fig. 3. The normalised free energy density / for (3 — 0.1. 



Because of the slow fall off at large s, the total free energy per unit length 
of the disclination line in infinite volume is logarithmically divergent. On the 
other hand, the integral of the difference of the free energies densities of 
the smooth and singular disclination lines is convergent, but it requires an 
estimate of the contribution of the core of the singular disclination line to 
the free energy. 

Let us estimate the free energy density of the singular disclination line. 
Outside the core the potential energy V(Q) has the ground state value ^ mm , 
and the elastic energy density given by the L\ term in formula (1) has the 
form 

8s 2 ' 

We assume that the core is formed when the total free energy density ap- 
proaches 0, because in the disordered phase Q = and then T vanishes. 
This means that the elastic energy at the boundary of the core is equal to 
l^minl- This condition gives the radius of the core 

R 

Furthermore, we assume that also inside the core the total free energy density 
vanishes. Then, the total free energy F c of the singular disclination line stored 
inside a concentric with it cylinder of unit height and of radius L£ is given 
by the integral 
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The factor 2tt comes from the integration over the angle 0. Note that F c 
includes the ground state energy ^ m i n . 

Now we can compare the free energy Fb of the smooth disclination line 
(inside the cylinder of radius L£ ), which is given by the formula 

rL 

F b = 2?r^ / dssF, 
Jo 

with F c . For the smooth disclination line we use the polynomial approxima- 
tion constructed in the previous Subsection. We have calculated numerically 
the difference Ft, — F c for (3 from 0.01 to 0.2 in steps of 0.01 taking L=3o. 
The results are plotted in Fig. 4. 



F b- F c 



-1 



-2 



0.01 0.05 0.1 0.15 0.2 

Fig. 4. Difference of the free energies of the smooth disclination line (Fb) and 
of the one with the singular core (F c ), in units ^f-arjoCo- 



We see that only for (3 < 0.12 the smooth disclination line has smaller 
total free energy. One can say that then the energetical cost of leaving 
the vacuum manifold is not too large. On the other hand, for larger (3 
it is energetically favourable to form the disordered core. Of course, our 
estimate of the "critical" value 0.12 of (3 depends on the assumptions about 
the formation and energy of the core. 

In Secton 3.2 we have noted that in the case of nematic MBBA (3 > 
0.174. Therefore, our calculations suggest that in that nematic material the 
n — \ disclination line has the singular core in the whole relevant range of 
temperature. 
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4 Discussion 



1. Let us recapitulate the main results of our work. We have found that 
there exists the special value 1/3 of the parameter (3 for which the smooth 
disclination line has particularly simple mathematical structure. We have 
obtained the approximate analytic formulas for the functions R, S. Finally, 
we have pointed out that for small (3 the smooth disclination line is ener- 
getically preferred, while for (3 large enough the singular core will appear. 
The critical value of (3 that we have obtained is approximately equal to 0.12. 
This result is based on several assumptions and, therefore, should rather be 
regarded as an estimate only. 

2. The approximate analytic description of the smooth disclination line pre- 
sented in subsection 3.3 has several advantages. First, at small and large 
s it approaches the exact solution by construction. The comparison with 
purely numerical solutions, see Figs. 1-2, shows that it is quite good also for 
intermediate values of s. Second, it has relatively simple form and it is easy 
to obtain. Moreover, it gives R, S in the full range [0, oo) of the independent 
variable s, as opposed to the numerical solution. Therefore, we think that 
in many cases the approximate analytic description can be quite a satisfac- 
tory substitute for the unknown exact analytic solution, as well as for the 
cumbersome purely numerical solution. Yet another argument for this comes 
from the fact that Landau-de Gennes model itself is also an approximation. 
Therefore, even the exact analytic or very precise numerical solutions of equa- 
tion (13) provide only approximate description of situation in a real nematic 
material. 

3. From physical viewpoint, perhaps the most interesting result of our paper 
is the suggestion that in the nematic MBBA the n = \ disclination lines 
have singular, disordered cores. However, let us remind our assumptions: 
L 2 = 0, Q = inside the core, the polynomial approximation for the smooth 
disclination line. In several papers, see, e.g., [19|, structure of defects in 
nematics confined in spherical droplets or in a capillary has been analysed 
in the framework of Landau-de Gennes theory. The defects correspond to 
hedgehogs or to n = 1 disclination lines. Various possibilities for the structure 
of the core have been found. We think that all these results are sufficiently 
interesting to motivate experimental investigations of the structure of the 
core. 
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